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Abstract 

We present two observations related to the application of linear (LFE) and nonlinear frac- 
tional equations (NFE). First, we give the comparison and estimates of the role of the fractional 
derivative term to the normal diffusion term in a LFE. The transition of the solution from nor- 
mal to anomalous transport is demonstrated and the dominant role of the power tails in the long 
time asymptotics is shown. Second, wave propagation or kinetics in a nonlinear media with frac- 
tal properties is considered. A corresponding fractional generalization of the Ginzburg-Landau 
and nonlinear Schrodinger equations is proposed. 



1 Introduction 



We call a fractional equation (FE) a differential equation that contains fractional derivatives or 
integrals. The awareness of the importance of this type of equation has grown continually in the 
last decade. Numerous applications have become apparent: wave propagation in a complex or 
porous media jjJ-H; random walks with a memory and flights j§-|6); kinetic theories of systems 
with chaotic 0-Q and pseudochaotic Jn| dynamics; and others (see j| and collections of papers 
in This new type of problem has increased rapidly in areas in which the fractal features of 

a process or the medium impose the necessity of using non-traditional tools in "regular" smooth 
physical equations. Exploitation of the fractional calculus for FE provides not only new types of 
mathematical constructions, but also new physical features of the described phenomenon 
While the linear FE has attracted fairly broad research activity, the study of nonlinear FE is at its 
very beginning |ll|-|Q. 

This paper consists of two parts on the FE properties. The first part is related to a linear FE 
(LFE) with the second (diffusional) derivative and the fractional derivative of order a. A goal of this 
part is to present a solution to the FE and compare roles of different parts of the FE. The second 
part is related to a nonlinear FE (NLFE) with cubic nonlinearity. We discuss some features of these 
equations, and possible applications. For a broad discussion of the fractal features of different media 
where waves propagate, see the elegant review ^]). In particular, a turbulent media of this type may 
be considered. We will show that such media lead to NLFEs in a natural way after simple change 
of the dispersion law. 



2 The interaction of normalized and anomalous transport 

The interest in and relevance of kinetic equations with fractional derivatives is a natural consequence 
of the realization of the importance of non-Gaussian properties of the statistics of many dynamical 
systems. There is already a substantial literature studying such equations in one or more space 
dimensions. In many cases of physical interest it would be reasonable to assume that both Gaussian 
and anomalous processes would play a role. One typically associates anomalous processes with 
algebraically decreasing tails of a probability distribution function (PDF) , while the bulk of the PDF 
is expected to be mostly Gaussian in chracter. In this note we explore the interaction of Gaussian 
and anomalous dynamics by means of a simple model kinetic equation in one space dimension for 
the PDF. More space dimensions would greatly complicate the analysis, but could readily be carried 
out. 
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We consider the kinetic equation with fractional derivatives 

dP(x,t) d 2 P d a P , 

-^F 1 = d^ + e W^' 1<a<2 (1) 

d a P 

where e is a constant, presumed small, but not necessarily so. The Riesz derivatives — — — is easily 

o|x'| Q 

d a P 

defined for our purposes in Fourier transformed space, so that the Fourier transform of — - — — is 

d\x'\ a 

— \k\ a P(k,t) where P(k,t) is the Fourier transform of P(x',t). We have not taken a fractional 
derivative with respect to time. We could do so, but this would be yet another complication in 
the analysis. We see that the Fourier transform of the right hand side of (Q) is —(k 2 + e\k\ a )P(,t). 
Thus, for large wavenumber and short wavelength the system exhibits normal, Gaussian transport, 
while for small wavenumber and large wavelength, the system is anomalous. The two phenomena 
are equal for wavenumber k? 

k T = e ^. (2) 

Thus ([!]) provides a simple model which mixes both normal and anomalous transport. 

It is convenient to change the independent variables from x',t' to x, t with the definitions 

x = e^x' (3) 

t = e&t' (4) 



so that 

dP(x,t) _ d 2 P d a P 



(5) 



dt dx 2 d\x\ a 

The coordinate (x' , t') are scaled to the Gaussian transport scale, while (x, t) are of a scale in which 
Gaussian and anomalous transport are comparable. We may easily represent a solution of (|J) which 
also satisfies 

P(x,0)=S(x) (6) 

and as a consequence 

POO 

P(x,t)dx = l, (7) 



as a Fourier integral: 

1 f°° 

P(x,t) = — dk exp[-t(k 2 + \k\ a ) - ikx}. (8) 

27r J-oo 

We compare (|^) with the solution of the equation with only anomalous transport 

dQ(x,t) d a Q{x,t) 



dt d\x\ a 

where Q also satisfies (@) and (^), so that 



(9) 



Q(x,t) = ^- / dkexp(-t\k\ a - ikx). (10) 

Z7T 
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The properties of Q(x,t) have been well-explored in the literature. We need such results for com- 
parison with the properties of P(x, i). 

We develop numerous series expansions of (0) and (JTc|) , some convergent and some asymptotic. 
When we know the behavior of P(x,t) and Q(x,t) for |a;| large, we shall also evaluate one spatial 
moment of P{x, t) and Q(x, t). With this information in hand we can then provide a description of 
the interaction of normal and anomalous transport. We start with (|^), which we rewrite as 



Re f°° 
P(x,t) = — / dkexp[- 

7T Jo 



t(k 2 + k a ) - ikx] 



We observe that 



where 



N 



-tk° 



y2(-m n —+ n N+1 (k,t), 



\R N +i{k,t)\ < t 



h.a(N+l) 
JV+1 K 



(JV+1)! 



(11) 



(12) 



(13) 



Since 



dk exp(— tk 2 — ikx) 



< 



t (N+l 



dk k^ N+1 h~ tk2 



(N + 1)1 Jo 

it is easy to conclude that the integral of the error tends to zero uniformly for all x and < ?/ < t < T 
as N — > oo, so that 



p(M) = -£(-i)^ / 



dk k an exp(-tk 2 - ikx), 



P(x,t) 



1 ^e-^ 784 ^, t n ^-^r{an+l) 

e~~ + Re — 5^(-l) n ^r^. —D- an -i 



2— r(n| 1) 



ix 



(14) 



(15) 



where D u (z) is the Weber function of order v. Again (15) converges uniformly for all x in < r\ < 
t < T. We see that the first term in ([l5]) exhibits Gaussian diffusion and also has unit integral. 
Hence, the remaining terms, the anomalous transport effects, have zero mean. For |x| large and n 
fixed we may employ the asymptotic expansion of the Weber function of large argument, or directly 
from (|l~5| ) we obtain the asymptotic expansion for large \x\: 



(-1) 



Tt + l 



n\\x\ 



1 4-n poo 

-Re(-l)"— / dk k an eM-tk 2 - ikx) 

b f<mn\ f ^ / t \ m T(an + 2m + l) 

— sin (— ) Eb 1 



M+l 



(16) 



The series ([16]) is clearly asymptotic and shows the power tail of the PDF, with the leading term 
from n = 1 being 

(t/\x\ a ) na 

!Alir(a + l)sm-. (17) 
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We may obtain another convergent expansion by use of the power series 

e- ikx = y2(-i) n ^-+R N+1 (k,x), (18) 
z — ' n! 

n=0 



where 



and thus 



where 



\r\ N+1 

i iiw+1 (A ; ,x)i<^ +i ^y T ^ I (19) 



P(x,t) = J2(~l) n x 2n C n (t), (20) 

n=0 



= TTTTi - / * k 2N exp-t(k 2 + k a ). (21) 
(2n)!7r 7 

Although C n (t) does ot appear to be one of the usual special functions it is easy to obtain an 
asymptotic expansion for t large and convergent expansion for t small. In particular for t small but 
bounded away from zero: 

C n (t) = ^ T ^V(-/^r (n + ™ + I) (22) 
" W (2n 2tt ^ m 12 2/ 1 ; 

v ; n=0 v 7 

and for i large 

We see from ( jlil) and (|l7| ) that for \x\ large P(x,t) ~ (const. )|x| _Q_1 . Thus, moments of order 
> a do not exist. We examine one, simple relevant moment 



M = \x\dx P(x,t) (24) 

J — oo 

and it follows easily that 

2 r°° r°° d 

M = — dxRe dk exp[-t(fc 2 + k a )] — sinfcz;. 
7T 7o Jo ' dk 

After integration by parts we obtain 

ni />oo />oo 

M = — / da; Im / dk(2k + ak^ 1 ) exp[-i(fc 2 + /fc Q )]e ifcr 
n Jo Jo 



(25) 



(26) 



If we move the path of x integration into the upper half-plane we find 



2t f°° 

M = — dk(2 + ak a - 2 )exp-t(k 2 + k a ) 
Jo 



With the same type of expansion as before we find an expansion convergent inO<rj<t<T 

'ma — 1 " 

-1)"'^ : 1' I 

m! \ 

m=0 



j.1/2 00 y.m(l-a/2) 

m = — V (-i) m+1 - . — r 

7r ^ — ' m! 



(27) 



(28) 
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The leading order term m — is just 2(t/n) 1 / 2 , exactly what one would obtain from a Gaussian. The 
remaining terms are the corrections from anomalous transport. For t large we have the asymptotic 
expansion 

M = fL—Y { - ir t- m ^-V [ ~ J \ a J (29) 
7r z — ' ml 

m=0 

We now turn to i) so that we may assess the significance of the preceding results. We start 
from ( JlOj ) and we find easily the convergent expansion 



1 00 TV 2n+l \ 



(30) 



while for |x| large there is the asymptotic expansion 



q — rr E(-ir(*/i*rr^_i) sin — (3 i) 

^Fl^ r(n + l) 2 

Finally, we define Mq as 

poo 

Mq= / Q(x,*)|x|<& (32) 



so that 



2 f°° [°° d 

Mn = — I a; da; Re / dk exp(— k a t)— sin kx. 
7T Jo Jo dk 



(33) 



We follow the same procedure as for M and we obtain 

M Q = -t^ a T(l - 1/a). (34) 



We are now prepared to compare results and we start with the simplest comparison, M and Mq . 
The result ([m]) for v is valid for all times, however M(t) is very different from (|34|) for t ~ 1. For t 
large (p9h shows 



7T 



iV-r ( i _ i) + ti-i/o-r fi + i ) + o(t 2 ~°) 



(35) 



Thus, in leading order and for 1 < a < 2, M ~ Mq. However, the difference between M and Mq is 
not small unless i 1 /" ^s> . For a not far from one, this relation holds, but as a approaches 2, one 
requires increasing larger values of t in order that t 1 /" dominate t 1_ i . It is clear that the expansion 



(29) fails at a — 2, and must then be very poor for small values of 2 — a. An examination of this 
moment indicates that for t sufficiently large the anomalous transport is the limiting form for the 
case with both Gaussian and anomalous transport, however there may be significant corrections. 
When we compare P(x, t) and Q(x, t) the situation is somewhat more complex. For t and x ~ 1, 



so that both Gaussian and normal transport may occur we may compare (|f 5[) , (Jl 6|) with (p0[),(|31|). 
When \x\ ~ f, P and Q are substantially different; however for t ~ 1, |sg| > 1, one may use (|f6|), 
and then the leading order terms in (|f6|) match M) although there are nontrivial corrections of order 
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t/x 2 . When t > 1 but |x| ~ 1, we may compare ( pi] ) and ( p3[ ) with (|30|). Again the leading order 
terms agree but there are nontrivial corrections of order t^d^. Finally, when both x and t are 
large, we must decide whether \x\ a /t is large or not. If |s| Q /i ~ 1 after a little effort we see that we 
may use (|2(]) and ( |23| ) for P(x,t), which again matches Q(x,t) in (|30|), but with corrections. Finally 
if |x| and t are both large, as is |a;| a /t, then we may conclude that ( fl4| ) and ([IB]) apply for P(x,t) 
which approximates Q(x,t) as given by (fjll). We see that broadly the anomalous transport finally 
dominates the Gaussian, or normal, transport, although it may take a long time. From another 
point of view, unless e is extremely small, the anomalous transport description is generally more 
relevant than the normal transport, although there may be significant corrections to the anomalous 
transport results. 

3 Fractional generalization of Ginzburg-Landau equation and 
nonlinear Schrodinger equation 

Let us recall the appearance of the nonlinear parabolic equation (see for example fL5| ). The simplest 
way is to consider a symmetric dispersion law lu — uj(k) for wave propagation in some media, and 
to represent the wave vector k in the form 

k = k + k = k + + k± (36) 

where ko is the unperturbed wave vector and subscripts (||, _L) are taken respectively to the direction 
ko. Considering k -C kg we have 

w{k) = wflko + k\) w w(fco) 

+ c(|k +k\- k ) rj w(fco) + ck„ + (37) 

2fc 



where c = duj/dk^. The expression ( |37| ) corresponds to the linear parabolic equation 

with respect to a field Z — Z(x,r,t), x along the kg, and the operator correspondence: 

d d d 

v = u(k) - u(ko) = igr, K \\ = ~ i ^ Kl - = r = (v> z ) ( 39 ) 

A generalization to a nonlinear case can be carried out in analogy with ( |37| ) through a nonlinear 
dispersion law depending on the wave amplitude: 

oj = oj(k, \Z\ 2 ) ^ uj{k,0) + g\Z\ 2 
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= W (|k + K|,0)+. 9 |Z| 2 (40) 

with some constant g = du>(k, \ Z\ 2 ) / d\Z\ 2 at \Z\ 2 = 0. In analogy with (|38|), the nonlinear parabolic 
equation takes the form 

dZ dZ c A , „ 

- i m= ic d^ + W AZ -^ z (41) 

This equation is also known as the nonlinear Schrodinger equation in which any of the coefficients 
may be complex. Indeed, for a traveling wave Z = Z[x — ct, r) we have 

-*# =i AZ -^i 2 < 42 > 

with £ = x — ct. 

Wave propagation in a media with fractal propeties can be easily generalized from the first step 
of writing dispersion law (^). Namely, one can replace ( |37| ) and (^) by the following: 

v = Lu{k, \Z\ 2 ) - u(ko, 0) = CK|| + ^(/^i)"/ 2 + g\Z\ 2 (43) 

with a fractional value of 1 < a < 2 and new constant c a . This replacement does not affect the 
nonlinear term which appears for K|| = k± — (see (|40|)). 

Using the connection between the differentiation operator and its Fourier transform 

(_A)«/2 <_> ( K a ja/a ( 44 ) 
we obtain the equation corresponding to ( ]43| ) : 

" ; f 4-5 ( -^ f -^ (45) 

or for travelling waves as in ( |42] ) 

iC f = + (46) 

Let us comment on the physical structure of (^) . The first term on the right-hand side is related to 
the wave propagation in a media with fractal properties. The fractal derivative may also appear as 
a result of ray chaos (l6) , or of a superdiffusive wave propagation (see also the discussion in [|| , Q 
and corresponding references). The second term on the right-hand side of (p5|),(p6[) corresponds to 
the wave interaction due to the nonlinear properties of the media. Thus, ( [46| ) describes fractional 
processes of self-focusing and related issues. 

We may consider a one-dimensional simplification of (^), i.e. 

"I = "H + ^' 2Z (47) 
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with some generalized constants c, a, g and then reduce (|47]) to the case of a propagating wave 
solution 

Z = Z(z-c£) = Z(ri). (48) 



Then (47) takes the form of fractional generalization of the Ginzburg-Landau equation (FGL): 

d a Z dZ „, 

a WF + ^ + 9Z =° (49) 

for real Z(rj). Now (^) differs from the fractional Burgers equation in the structure of 

the nonlinear term. Nevertheless, an analysis similar to |]l3| , ||l4]] may be performed to obtain some 
estimates on the solution. 

It is well known that the nonlinear term in the equations of ( ff2| ) type leads to a steepening of the 
solution and its singularity. The steepening process may be stopped by a diffusional or dispersional 
term, i.e. by a higher derivative term. A similar phenomenon may appear for the fractional nonlinear 
equations (f47j) , (ff9|) . It has been shown in |l3| that for the fractional Burgers equation there exists 
a critical value a c such that solution is regular for all time if a > a c . 

Next, we consider the symmetric solution Z{rj) = Z(—rf) of ( ]48| ) and assume Z € IR 2 . Let us 
multiply ( |49| ) by Z and integrate it under the assumption that Z is sufficiently small at infinity. We 
find 

/>oo 

£= / {aZZ (a ^ +gZ A )dri = (50) 



where we use the notation 

£t- = zia) ( 51 ) 

d\v\ a 



The expression (50) gives a conservation law showing the relative importance of the nonlinear and 
dispersion terms. Considering Z^ ~ \jrp we conclude that the term with derivative prevails for 
small rj if 

a > 2 7 . (52) 



For 2 > a > 1 the condition (52) shows that the solution is square integrable at rj — > 0. 

Finally we observe that (|4^) may describe a new type of fractional dynamic motion. This topic 
will be considered elsewhere. 



4 Conclusion 

In this brief paper we have made two remarks related to application of fractional equations in physics. 
The first remark is related to the competition between normal diffusion and diffusion induced by 
fractional derivatives, as occur in fractional kinetic theories. It is shown that for large times the 
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fractional derivative term dominates the solution and leads to power type tails. The second remark is 
related to a new class of equations in which fractional derivative terms are responsible for fractional 
dispersion. In this case the asymptotics should be defined by a competition between the fractional 
dispersion and nonlinear terms. We discuss the origin of the fractional Ginzburg-Landau equation 
and fractional nonlinear Schrodingcr equation. 
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